Spin-orbit interactions in black-hole binaries 
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We perform numerical simulations of black-hole binaries to study the exchange of spin and orbital 
angular momentum during the last, highly nonlinear, stages of the coalescence process. To calculate 
the transfer of angular momentum from orbital to spin, we start with two quasi-circular configu- 
rations, one with initially non-spinning black holes, the other with corotating black holes. In both 
cases the binaries complete almost two orbits before merging. We find that, during these last orbits, 
the specific spin (a/m) of each horizon increases by only 0.012 for the initially non-spinning configu- 
ration, and by only 0.006 for the initially corotating configuration. By contrast, the corotation value 
for the specific spin should increase from 0.1 at the initial proper separation of 10M to 0.33 when 
the proper separation is 5M. Thus the spin-orbit coupling is far too weak to tidally lock the binary 
to a corotating state during the late-inspiral phase. We also study the converse transfer from spin 
into orbital motion. In this case, we start the simulations with parallel, highly-spinning non-boosted 
black holes. As the collision proceeds, the system acquires a non-head-on orbital motion, due to 
spin-orbit coupling, that leads to the radiation of angular momentum. We are able to accurately 
measure the energy and angular momentum losses and model their dependence on the initial spins. 

PACS numbers: 04.25.Dm, 04.25. Nx, 04.30.Db, 04.70.Bw 



I. INTRODUCTION 



With the advent of new methods to evolve black-hole 
binaries for several orbits 0, d, d, 0, H, H, @] we can 
now pose questions of physical interest about the details 
of their interactions and the generation of gravitational 
radiation, and find answers that allow us to better under- 
stand the astrophysical properties of black-hole binaries. 

The coalescence of unequal-mass black-hole binaries 
with randomly oriented individual spins is expected to 
be one of the strongest astrophysical sources of gravita- 
tional radiation. The first numerical simulations of non- 
spinning unequal-mass black- hole binaries, along with 
the calculation of the merger kicks, were reported in 
Refs. USES. 

In a preceding paper [Tl| we performed the first fully- 
nonlinear simulations of highly spinning black- hole bina- 
ries and showed that the spin-orbit coupling leads to sig- 
nificant changes in the orbital dynamics, as well as grav- 
itation radiation emitted, in the final stages of a black- 
hole binary. In this paper we study in more detail the 
mechanisms that transfer orbital angular momentum to 
spin and vice versa. 

It is of astrophysical interest to know if tidal effects 
could lock the spin of black-hole binaries to a corotating 
state, i.e. a state where the spins of the individual hori- 
zons are aligned with orbital angular momentum, and 
the spin magnitudes are such that the horizon frequency 
£Ih (Oij is the angular speed of locally-non-rotating ob- 
servers as they pass through the horizon, as seen by sta- 
tionary observers at infinity) is nearly equal to the orbital 
frequency Qo- For a Kerr hole with mass m and specific 



spin a/m the horizon frequency takes the value 



Ira 



a/m 



1 + yj\ - (a/m) 2 



(1) 



See Ref. [12j for a first post-Newtonian (1PN) correction 
to the corotation condition f2<p = Note that in this 
paper we consider both the spin (or more accurately, the 
spin angular momentum) S of a black hole, which has 
units of M 2 , and the specific-spin a/m (m is the horizon 
mass), which is dimensionless. The two are related by 
S = am. We will also make use of the following relation- 
ships between the spin, specific spin, horizon frequency, 
and irreducible mass (m IR ): 



S = 



y/i - i(m 1R n H ) 2 ' 



a/m = 4m 1R tt H \/l - A(m m tt H ) 2 , 



(2) 
(3) 



where 167rrn 2 R is the horizon surface area. 

There are two competing scales in this problem: one 
is the time for tidal effects to act to spin up (or down) 
the black holes to the orbital frequency and the other is 
the rate of change of the orbital frequency due to gravi- 
tational radiation. 

In the early stages of the orbital motion of black-hole 
binaries, tidal effects are small and can be described by 
Newtonian theory; these effects decay strongly with the 
binary separation r (as /i/r 3 , with \x the reduced mass of 
the system). Tidal effects will be stronger as the holes 
get closer, but this will also increase the orbital frequency 
(roughly like f^o ~ \J /V 7 " 3 )- In order to be correctly 
quantified, these effects need to be evaluated dynami- 
cally with a fully non-linear numerical simulation of the 
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final orbital stage of black-hole binaries. Despite this, 
several computations of quasi-circular sequences of ini- 
tial data for black-hole binaries have assumed corotat- 
ing spins from larger separations down to the estimated 
location of the innermost stable circular orbit (ISCO) 
E3) EH • Corotation has a notable effect on the lo- 
cation of the ISCO in all these calculations; fortunately, 
recent numerical evolutions of equal-mass black-hole bi- 
naries 0, [H, [f| do not find any particular physical effect 
associated with the location of the ISCO. 

Bildsten and Cutler [l5[ (see also studied neutron 
star (NS) and black hole (BH) systems. They conclude 
that, for the NS-NS binary the viscosity is too low to 
tidally synchronize the stars as they inspiral, and that 
there is simply not enough time during the gravitational 
radiation driven inspiral of NS-BH binaries to tidally lock 
the spin of the neutron star to the orbital period. Addi- 
tionally, when the NS-BH binaries are close enough for 
the tidal effects to be strong, the neutron stars are dis- 
rupted in a few orbits. Since these arguments do not 
strictly apply to black-hole binaries, it is interesting to 
explore these tidal interactions dynamically with fully 
non-linear numerical techniques. This spin-up of NS-NS 
binaries has also been analyzed in Ref. [17] . 

We will study in SeaQHthe tidal spin-up of two equal- 
mass, initially non-spinning, black holes starting from a 
quasi-circular orbit (as determined by 3PN trajectories) 
with nearly two orbits left before merger, as well as a 
similar configuration with initially corotating black holes. 
These equal-mass configurations involving the very last 
stages of the inspiral should nearly maximize the tidal 
effects. 

While the transfer of orbital angular momentum to the 
spin of the members of a binary system is of astrophysical 
relevance, we can also study the conversion of spin into 
orbital angular momentum and verify the importance of 
its dynamical effects. To this end, we evolve binaries con- 
sisting of non-boosted, equal- mass black holes, with the 
black holes aligned along the y-axis and the individual 
spins aligned along the +z-axis, and measure the radi- 
ated energy and angular momentum. This simple set-up 
clearly displays the transfer of spin to orbital motion. 
For non-spinning holes this set-up produces purely head- 
on collisions, leading to the radiation of energy but not 
angular momentum. For spinning holes, the spin-orbit 
coupling causes the holes to acquire non-head-on orbital 
motion, leading to the radiation of angular momentum. 
Thus, an accurate measure of the angular momentum ra- 
diated would allow us to estimate the correct magnitude 
of this effect. 

The leading post-Newtonian spin-orbit and spin-spin 
corrections to the acceleration of two equal-mass black 
holes (each of mass m), with spins Si and S2 perpendic- 
ular to the motion, in a nearly head-on collision are: 

16 - 
aso = — 3 fa x (Si + 5" 2 )), (4) 

6 - 

ass = j{Si ■ S 2 ) f, (5) 



where r is the binary separation, v r is the relative ve- 
locity, and SO and SS denote spin-orbit and spin-spin 
coupling terms respectively. While the spin-spin acceler- 
ation adds an attractive term along the radial direction 
(for parallel spins; repulsive for anti-parallel), the spin- 
orbit acceleration is perpendicular to both the mostly- 
radial (along r) head-on velocity and the spin, and hence 
leads to a non- vanishing orbital component of the binary 
motion. 

In Sec. IIVI we will study this nonlinear coupling us- 
ing fully non-linear numerical evolutions and show that 
our method is able to extract meaningful results for this 
numerically demanding relatively small effect. 



II. TECHNIQUES 

We use the Brandt-Briigmann puncture ap- 
proach [ll| along with the TwoPunctures [l9 ] 
and BAM_Elliptic [2(j thorns to compute initial 
data. In this approach the 3-metric on the initial 
slice has the form 7^ = (ipBL + u) 4 °~ab, where ipBL 
is the Brill- Lindquist conformal factor, 5 a t, is the 
Euclidean metric, and u is (at least) C 2 on the punc- 
tures. The Brill-Lindquist conformal factor is given by 
ijjBL = 1 + Y^j—i Ttii/(2ri), where n is the total number 
of 'punctures', m.j is the mass parameter of puncture i 
(rrii is not the horizon mass associated with puncture 
i), and r, is the coordinate distance to puncture i. In 
all cases below, we evolve data containing only two 
punctures with equal puncture mass parameters, and we 
denote this puncture mass parameter by m p . We evolve 
these black-hole-binary data sets using the LazEv [21| 
implementation of the moving puncture approach [2, S] . 
In our version of the moving puncture approach [2| we 
replace the BSSN [2^, HH, [2J] conformal exponent (j>, 
which has logarithmic singularities at the punctures, 
with the initially C 4 field \ = exp(— 40). This new 
variable, along with the other BSSN variables, will 
remain finite provided that one uses a suitable choice for 
the gauge. An alternative approach uses standard finite 
differencing of (f> ||. Note that both approaches have 
been used successfully by several other groups d, H^, . 

We obtain accurate, convergent waveforms and horizon 
parameters by evolving this system in conjunction with 
a modified 1+log lapse, a modified Gamma-driver shift 
condition [2|,[27|, and an initial lapse a ~ ^bl- The l a P se 
and shift are evolved with (dt — f3 l di)a = —2aK, dt/3 a = 
B a , and dtB a = i/AdtT a — nB a . These gauge conditions 
require careful treatment of \ near the puncture in order 
for the system to remain stable @, 0). In Ref. [28| it 
was shown that this choice of gauge leads to a strongly 
hyperbolic evolution system provided that the shift does 
not become too large. For our version of the moving 
puncture approach, we find that the product aA ll dj(f> 
initially has to be C 4 on the puncture. In the spinning 
case, is 0(r 3 ) on the puncture, thus requiring that 
a oc r 3 to maintain differentiability. We therefore choose 
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an initial lapse a(t = 0) = 2/(1 + ip% L ) which is 0(r 4 ) 
and C on the puncture and reproduces the isotropic 
Schwarzschild lapse at large distances from the horizons. 
The initial values of fj % and B % are set to zero. 

Hannam et. al. [261 ] examine the smoothness of the 
evolved fields at late times at the puncture. They find 
that, in the case of Schwarzschild, x transitions from an 
initially C 4 field to a C 2 field at late times. Although 
we require that the fields are initially C 4 , this late-time 
drop in smoothness does not appear to leak out of the 
horizon (which is consistent with the analysis in 26]). 

We use a 'multiple transition' fisheye transforma- 
tion [3] to push the boundaries to lOOAf, while main- 
taining a resolution of up to M / 40 in the central region, 
for the head-on runs; for orbital runs we push the bound- 
aries to 215M, with a central resolution of up to M/31.5. 

We use a locally modified version of the AHFinderDi- 
rect thorn [1^] to find the individual apparent horizon, 
and we measure the angular momentum of the horizons 
using our implementation of the isolated horizon algo- 
rithm detailed in [3(J. (The isolated horizon formalism 
used here is extended to the dynamical case in |3l|.) In 
the isolated horizon approach, the horizon spin is given 
by 

S = ^-<f (p a R b K ab )d 2 V. (6) 

^ J AH 

where ip a is an approximate Killing vector on the horizon, 
K a i, is the extrinsic curvature, d 2 V is the natural area 
element intrinsic to the horizon, and R a is the outward 
pointing unit vector normal to the horizon on the 3D 
slice. (See appendix [A] for a comparison of the isolated 
horizon technique for finding the horizon spin to the more 
common horizon circumference technique.) We confirm 
that our version of the algorithm produces the correct 
spin for various coordinate distortions of Kerr, as well 
as for Bowen-York binary initial data (both with and 
without linear momentum). 

We measure the proper distance between the horizons 
by shooting geodesies from the origin to a horizon, and 
then finding the geodesic with the minimum proper dis- 
tance as described in Refs. [12, [HI (the proper distance 
between the two holes is then double this distance) . This 
technique will work in cases where there is (at least) Pi- 
symmetry on the equatorial plane and reflection symme- 
try across the plane. 

III. TIDAL SPIN-UP OF BLACK-HOLE 
BINARIES 

The spin of the individual horizons in a black-hole bi- 
nary will become tidally locked to the orbital frequency 
if the non-linear spin-orbit coupling is sufficiently strong 
that the timescale for tidally induced spin-up (or down) 
is smaller than the timescale for change in the orbital 
period. In order to obtain an accurate measure of the 
strength of this coupling we examine the cases of initially 



TABLE I: Initial data for quasi-circular, equal-mass black- 
hole binaries. The binaries have an ADM mass of (1.0000 ± 
0.0005)M, with orbital frequency MQ. fixed to 0.0500, total 
angular momentum (orbital + spin) given by (0,0, J), and 
proper horizon separation I. The punctures are located at 
(0, ±Y, 0), with mass parameter m p , momentum (+P, 0, 0), 
spin angular momentum (0, 0, S), and specific spin a/m (m is 
the horizon mass). 

Name S/M 2 Y/M P/M J/M 2 a/m m p /M l/M 
SO O000 3.280 0.1336 0.876 0.0000 0.4848 10.01 
SC 0.025757 3.2534 0.1330 0.917 0.1001 0.4831 9.93 



non-spinning and initially corotating close binaries, de- 
noted by SO and SC respectively, where the binary sep- 
aration is small enough that the spin-orbit coupling is 
large, but large enough that the binaries complete at least 
~ 1.75 orbits before merging. The initial data parame- 
ters for these configuration are provided in Table U and 
have been obtained by finding quasi-circular orbits with 
the 3PN equations of motion. We originally studied the 
SO non-spinning configuration in Ref. [ill ] . While in [TTI ] 
we were interested in the spin of the remnant horizon 
and plunge times as a function of the initial spin of the 
individual horizons, here we are interested only in the 
pre-merger spin-up of the individual horizons. It is then 
crucial to have an appropriate measure of the individ- 
ual spin of the holes. We employ the isolated horizon 
spin algorithm detailed in [30] (see Eq. (BJ). Strictly 
speaking, isolated horizon techniques are only valid for 
non-interacting horizons. However, they do provide an 
excellent approximation to the horizon spin even in the 
case of relatively close binaries. We confirmed that our 
implementation gives accurate evaluations of the horizon 
spins for closely separated Bowen-York binary data. For 
this test we constructed a sequence of initial data using 
the TwoPunctures thorn [lj| with fixed orbital angu- 
lar momentum L = 0.876M 2 (the value for the SO case 
above) and ADM mass of (1.0000±0.0005)M. We set the 
spin of the individual holes (aligned along the +z-axis) 
to S = OM 2 , S = 0.001M 2 , and S = 0.0001M 2 , and 
measured the spins of the horizons with our implemen- 
tation of the isolated horizon algorithm. The results of 
these calculations are summarized in Table [TTJ From the 
table we see that spins as low as 10 _4 M 2 can be mea- 
sured with an accuracy of better than 20% even when the 
holes are separated by a proper distance of / ~ 3M, and 
can be measured with an accuracy of 1% or better when 
the holes are separated by I > 4.6M. Interestingly, ac- 
curate measurement of this spin can be obtain even with 
coordinate separations as small as 0.5M. 

We evolved the SO and SC configurations with 3 dif- 
ferent resolutions, h = Af/22.5, h = M/27, and h = 
M/31.5, where h is the central grid-spacing, and calcu- 
lated the horizon spins after every 1M of evolution. The 
results of measuring the spin of the individual holes dur- 
ing the orbital motion until the formation of a common 
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TABLE II: The isolated horizon spin of a sequence of 
Bowen-York data with given spins and ADM mass (1.0000 ± 
0.0005)M. Here (±x,0,0) are the coordinate locations of the 
punctures, S is the measured spin (the Bowen-York spin pa- 
rameter is given in parentheses), and I is the proper distance 
between horizons. In all cases the linear momentum is given 
by P = (0, ±0.8764Af 2 /(2a;),0) (to reproduce the SO orbital 
angular momentum), and the spins are aligned with the z- 



axis. 


x/M 


S/M 2 (KT 3 ) 


S/M 2 (10" 4 ) 


S/M 2 {0) 


l/M 


3.280 


1.000 x 10~ 3 


1.003 x 10" 4 


3.0 x 10" ' 


10.0 


3.0 


1.011 x 10~ 3 


1.011 x 10~ 4 


1.1 x 10~ 6 


9.4 


2.0 


9.998 x 10" 4 


9.97 x 10" 5 


3.4 x 10~ 7 
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1.0 
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0.5 
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TABLE III: Merger time, T C ah/M, for the SO and SC config- 
urations versus resolution. Horizon searches were performed 
every 0.3M and 0.2M for the SO and SC configuration respec- 
tively. 



resolution 


SO 


SC 


M/22.5 


160.7 ± 0.3 


168.6 ±0.2 


M/27 


166.0 ± 0.3 


174.2 ±0.2 


M/31.5 


168.3 ± 0.3 


176.6 ±0.2 


extrapolation 


172 ±2 


179 ±2 



apparent horizon (CAH) (See Table IIII[) for the SO con- 
figuration are displayed in Fig. [TJ while the results for 
the SC configuration are displayed in Fig. [2l The up- 
per panel in both figures shows the specific spin of the 
individual holes for the three resolutions. We observe 
that in both cases the specific-spin increase during the 
evolution is small (0.012 for SO and 0.006 for SC). The 
lower panels show essentially third-order convergence as 
expected from our use of third-order-accurate interpola- 
tors with AHFinderDirect. The late-time error in the 
SO specific spin for the lowest resolution run appears to 
show a very large error compared to the actual value of 
the specific spin. This 'error' is actually a phase error 
resulting from the lower resolution runs merging sooner. 
However, in a specific spin versus proper separation plot, 
the curves from all three resolutions lie on top of each 
other (see Fig. [3]) at late-times (i.e. at small proper dis- 
tance). This dichotomy, between the late-time conver- 
gence break-down of a/m versus t and the excellent late- 
time agreement between resolutions in a/m versus I, is 
similar to the resultspreviously seen in the waveforms of 
inspiraling binaries 0, [111, |34j. In that case it was ob- 
served that iplit) lost convergence due to large phase er- 
rors induced by the lower resolution runs merging sooner, 
while the phase-corrected ip4 showed excellent agreement 
between resolutions. In this case the 'phase error' is in 
l(t), and by plotting a/m versus I we effectively removed 
the phase error. 

The values that we obtained for the spin-up of the 



FIG. 1: The measured instantaneous specific spin of the indi- 
vidual horizons, for the SO configuration, starting from van- 
ishing spin at t = and up to the merger time. The top 
panel shows a/m versus time for three resolutions with grid- 
spacings h = M/22.5, h = M/27, and h = M/31.5, while 
the lower panel shows the differences in a/m between the low 
and medium resolutions and the medium and high resolutions 
(the latter rescaled by 1.9662 to demonstrate third-order con- 
vergence). The differences between the low and medium reso- 
lutions and medium and high resolutions become large at late 
times due to the lower resolution runs merging sooner. The 
curves in the top panel have been cut off at the approximate 
merger value of a/m (which is independent of resolution). 



binary holes are much smaller than those expected for 
a corotation state [43j |. In order to visualize this we 
have displayed in Fig. U the measured spin during the 
orbital motion, this time against the proper separation 
of the apparent horizons. The SO configuration forms 
a common apparent horizon (CAH) at a proper separa- 
tion of Icah — 1.75M, while the SC configuration forms 
a common apparent horizon at a proper separation of 
Icah = 1.69M. In order to have a measure of what the 
corresponding corotation spin is, we considered the fam- 
ily of corotating binaries in quasi-equilibrium of Ref. [lij 
(we have used Eq. ([3]) to convert the and rrij pro- 
vided in [lj| into a/m). From the plot we see that the 
SC is initially corotating, but the specific spin remains 
flat for most of the run. Thus the spin-orbit coupling 
is far too small to keep this system tidally locked. It is 
interesting to note that the spin-up for the SO configura- 
tion is slightly larger than the spin-up for the SC config- 
uration. Both configuration show essentially no spin-up 
above proper separations of ~ AM . 

The relatively low 'spin up' that we measure during 
the last couple of orbits prior to the merger indicates that 
the time-scale for tidal locking is much larger than the 
time-scale of orbital decay due to gravitational radiation 
emission. We also note that in the leading relevant post- 
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FIG. 2: The measured instantaneous specific spin of the in- 
dividual horizons, for the SC configuration, starting from 
a/m = 0.1 at t = and up to the merger time. The top 
panel shows a/m versus time for three resolutions with grid- 
spacings h = M/22.5, h = M/27, and h = M/31.5, while 
the lower panel shows the differences in a/m between the low 
and medium resolutions and the medium and high resolutions 
(the latter rescaled by 1.9662 to demonstrate third-order con- 
vergence). The differences between the low and medium res- 
olutions and the medium and high resolutions become large 
at late times due to the lower resolution runs merging sooner. 
The curves in the top panel have been cut off at the approx- 
imate merger value of a/m (which is independent of resolu- 
tion) . 
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FIG. 3: The specific spin a/m for the SC configuration versus 
proper binary separation / for three resolutions with grid- 
spacings h = M/22.5, h = M/27, and h = M/31.5. Note 
that the large late-time phase errors seen in a/m versus time 
are not present in this plot. The vertical line shows the proper 
distance at merger. 
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FIG. 4: The spin-up of the horizons versus proper binary 
separation. The continuous line shows a/m for the SO config- 
uration, the dotted line shows a/m for the SC configuration 
(initially corotating) , and the dot-dash line shows the corota- 
tion value for the spin. Note that the SO configuration shows 
a slightly larger spin-up and that both spin-ups are much 
smaller than what is required to tidally lock the binary. 



Newtonian order, the change of the spin, averaged over 
a circular orbit, which is given by [35] ] 



±[(L N xSi)(4 + 3^) + 5 J x5 i 



-3(L N ■ §j)L N x (j ? i), 



(7) 



vanishes for all configurations with the spins aligned with 
the orbital angular momentum. In Eq. ([7]) Ln is the 
(Newtonian) orbital angular momentum, Si is the spin 
of black hole i, rrii is its mass, and j denotes the other 
black hole. Will [36j also points out that the radiation 
reaction interaction does not modify the spin up to 3.5PN 
order. First effects of the flattening of the bodies should 
actually be seen at 4.5PN order. 



IV. SPIN-ORBIT COUPLINGS IN 
BLACK-HOLE BINARIES 

We continue our study of the relevance of spin-orbit 
coupling to the orbital dynamics of black-hole binaries 
by choosing initially head-on (vanishing transverse linear 
momentum) configurations that single out the converse 
transfer of spin into orbital angular momentum. Purely 
head-on collisions do not radiate angular momentum to 
infinity. Thus, in these configurations, the radiated angu- 
lar momentum is an indicator of the extent to which these 
systems deviate from head-on motion (i.e. the extent to 
which they acquire transverse motion). The energy and 
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TABLE IV: Initial data for non-boosted black-hole binaries. 
The binaries have an ADM mass of (1.0000 ± 0.00005)M, 
with total angular momentum (0,0, J), and proper horizon 
separation I. The punctures are located at (0,±Y,0), with 
mass parameter m p , spins angular momentum (0,0,5*), and 
specific spin a/m (here m is the horizon mass). 



Name 


S/M 2 Y/M 


P/M J/M 2 


a/m 


rrip/M 


l/M 


HS++ 


0.200 4.000 


0.000 0.4 


0.7581 


0.3444 


12.24 


MS++ 


0.100 4.2422 


0.000 0.2 


0.3775 


0.469755 


12.24 


0S++ 


0.000 4.2924 


0.000 0.0 


0.0000 


0.5000 


12.24 



angular momentum radiated are functions of the spins of 
the two horizons, the initial inward momentum of each 
hole, and the proper separation. To simplify the anal- 
ysis, we choose head-on binaries with zero initial linear 
momentum and with a fixed initial proper separation of 
12.24M. All spins are aligned along the z-axis, with the 
binaries initially separated along the y-axis. We report 
the initial data parameters in Table [TV] 

We measure the radiated energy and angular momen- 
tum by solving for the final hole remnant horizon mass 
and spin, and then subtracting these quantities from the 
ADM mass and angular momentum (as determined from 
the initial data). To avoid confusion with the initial 
black-hole mass and spins, we denote the mass and spin 
of the final remnant horizon with Mh and Jh respec- 
tively. The radiated mass and angular momentum then 
are: 



E ra d = MADM ~ Mjj, 
J rad — J ADM — Jh, 



(8) 
(9) 



where in these head-on cases J adm is the sum of the 
spins of the two initial black holes. For notational con- 
venience we introduce the symbol jo defined by 



Jo 



Jadm 

ADM 



(10) 



In previous papers [3, [TTJ] we found that the horizon 
provided a very accurate measurement of these quan- 
tities. Additionally, when using isolated horizon tech- 
niques, we can obtain the horizon parameters accurately 
within 30M after the merger. In all cases the radiated 
energy and angular momentum are very small. In order 
to accurately obtain these numbers we need to evolve 
all configurations at three different resolutions (see Ta- 
ble [V] for the grid-spacings) and then use a Richardson 
extrapolation of the radiated energy and angular momen- 
tum. The measured convergence rates for the radiated 
angular momentum are 3.9 and 3.7 for the HSH — h and 



MS++ configurations respectively, while the measured 
convergence rates for the radiated energy are 4.0, 3.7, 
and 3.9 for the HS++, MS++, and 0S++ configurations 
respectively. The lower convergence rates for the MS+H- 
configuration may simply indicate only that higher res- 
olution is needed. Alternatively it is possible that small 
third-order errors (possibly from the interpolators) are 
present in all configurations, but are only observed in 
the MS++ configuration due to its finer effective resolu- 
tion. (The effective grid-spacing is h e ff — h/m p , where 
h is the actual grid -spacing and m p is the puncture mass 
parameter (see for details)). For the Richardson ex- 
trapolation we use both the forms f(h) = /oo + e h 4 +g h 5 
and f(h) = foo + dh 3 + e /i 4 , and, using results from the 
three resolutions, solve for f^. We take the difference 
between these two extrapolation to be indicative of the 
error in our extrapolated numbers. In Figures [5] and [6] 
we show the HSH — h final remnant black-hole horizon mass 
and spin versus time for three resolutions, as well as the 
two extrapolations of these quantities and convergence 
plots. It is evident from the very close agreement be- 
tween the two extrapolations that we obtain highly ac- 
curate measurements of both the horizon mass and spin, 
and consequently, the radiated mass and angular momen- 
tum. 

In Table [V] we summarize the results for these con- 
figurations. Note that the errors reported in the table 
for a given resolution are of the uncertainty in measure- 
ment of Jh and Mjj at that resolution (i.e. a measure 
of the flatness of the late-time plateaus of Jh and Mr). 
Thus the error interval reported for e.g. J radiated at a 
resolution of M/25 will not necessarily contain the final 
extrapolated value. On the other hand, the error inter- 
vals for the extrapolated values are our best estimate for 
the true value of each quantity. A fit of the Richardson- 
extrapolated values of the radiated angular momentum 
to the initial spins of the individual black holes has the 
form: 



J, 



rad 



AP 



(4.8 ±0.1) x 10~ 3 j + 

(4.0 ± 0.9) x 10- 3 j 3 + O (jo 5 ) , (11) 



where j = 2S/M 2 ADM (see Table HB 



The linear term 
is the dominant one for these configurations (the HS+H- 
configuration radiates ~ 20% more momentum than ex- 
pected from the linear behavior alone). Equation (fTT|) 
contains only odd powers of jo since the magnitude of 
Jrad must be independent of jo (i.e. the physics would 
not change if the two spins were aligned along the —z- 
axis), and the sign of J ra d must match the sign of jo. 



An analysis of the radiated energy versus spin is com- 



plicated by the fact that the binary radiates energy in 
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TABLE V: The radiated mass and angular momentum as well as the merger time (Tcah) and final specific spin for the head-on 
configurations. The reported error intervals in Tcah are a function of how often a horizon search was performed. In all cases 
Tcah does not appear to vary with resolution. 



Case 


resolution 


E ra d/ MADM 


Jradj J ADM 


Tcah/M 


a/M H 


HS++ 


M/25 


(0.16 ±0.02)% 


(2.37 ±0.02)% 


36.7 ±0.2 


0.3918 ± 0.002 




M/30 


(0.136 ±0.009)% 


(1.44 ±0.01)% 


36.6 ±0.2 


0.39532 ± 0.00008 




M/40 


(0.123 ±0.004)% 


(0.832 ± 0.003)% 


36.7 ±0.2 


0.39765 ± 0.00003 




Avf ~r o t~\ f \ mt"i r^n 
(JXLI aJJUlclllUll 


(ft 1 1 8 4- ft ftft9\°/r, 


(ft ^ztR 4- ft nn9^% 


qc 7 _i_ ri n 
OU. 1 HI U.Z 


n iqr7r 4- n nnnn9 


MS++ 


M/25 


(0.065 ± 0.004)% 


(1.038 ± 0.008)% 


39.5 ±0.2 


0.19817 ±0.00003 




M/30 


(0.063 ±0.003)% 


(0.7698 ± 0.0040)% 


39.6 ±0.2 


0.19868 ±0.00003 




M/35 


(0.062 ±0.002)% 


(0.6567 ± 0.0069)% 


39.5 ±0.2 


0.19891 ±0.00003 




extrapolation 


(0.060 ±0.002)% 


(0.4986 ± 0.0086)% 


39.5 ±0.2 


0.19920 ±0.00007 


0S++ 


M/22.5 


0.0616 ± 0.0032 





40.4 ±0.2 







M/25 


(0.0595 ± 0.0025)% 





40.4 ±0.2 







M/30 


(0.057 ± 0.002)% 





40.4 ±0.2 







extrapolation 


(0.054 ± 0.002)% 





40.4 ±0.2 






the zero-spin configuration as well. After subtracting the 
zero-spin radiated energy we find that the HS±± config- 
uration radiates 0.064 ± 0.003M more energy than the 
zero-spin configuration, while the MS-I — h configuration 
only radiates 0.006±0.002M more. Thus the HS±± con- 
figuration radiates 12.2±4.6 times as much as the MS±± 
configuration (after subtracting off the zero-spin radiated 
energy). From these results it appears that the radiated 
energy scales with j$ or jg. A close-limit analysis (See 
Eq. |H])) also indicates that there is a small j'q contribu- 
tion as well. Although not obvious in the total radiated 
energy, this quadratic dependence on J can be seen in the 
imaginary part of the {I = 2,m = ±2) components of the 
radiated energy E(t) (i.e. the amount of radiation emit- 
ted up to time t), where E{t) = J2e, m ( R im(t) + hm{t)), 

1 r* 

Rlm{t) = — / Re[N im (r)] 2 dT, (12) 

Itm(f) = t / Im[N tm (T)} 2 dT, (13) 
4tt Jo 

Ne m (t) = (r - 2M)/2 J Q in,i m {r) dr, and r/> 4 is calcu- 
lated using the quasi-Kinnersley tetrad method [37] . Fig- 
ure [7] shows Re m (t) and Ii m {t) for the I = 2 components 
of the waveform for the HS±± and MS±± configurations. 
Note that I 2 ±2 (which vanishes in the zero-spin case) is 
4.3 times larger in the HS±± configuration than in the 
MS-I — h configuration. This is consistent with a quadratic 
dependence on the spin plus smaller higher-order depen- 
dencies. 

Finally, if we fit the radiated energy to the polynomial 
dependence E + E 2 j 2 ± we obtain 

= (0.052 ± 0.02) ± (0.065 ± 0.033)j 2 

±(2.08±0.21)j 4 ±O(j 6 ). (14) 

Note that the coefficient of the quadratic term is two 
orders of magnitude smaller than the coefficient of the 
quartic term (and thus was not apparent in our initial 



analysis above). There are only even terms in Eq. (|14|) 
because the E ra ^ cannot depend on the sign of j Q (i.e. 
E ra d must be the same for configurations with the spins 
aligned on the ±z-axis and with the spins aligned on the 
— z-axis). 

The Bowen-York initial data family contains 'spurious' 
radiation that does not represent any physical history of 
the binary. In addition, the spinning case does not have 
the Kerr limit when the two black holes are far apart. 
(For a data set with the right Kerr limit see [H, |39|.) 
While this extra radiation is of relatively small magnitude 
in orbiting black- hole binaries (See ) , it is relevant 
in the highly spinning head-on collisions considered here. 
This can be observed in Fig. [7) The first plateau in the 
radiated energy (t = 40M — 70M) for the observer loca- 
tion at r = 35M is due to the radiated 'spurious' energy. 
This 'spurious' energy contributes up to 25% of the total 
energy radiated of the HS-I — I- configuration, but strongly 
decreases for the MS±± and 0S±± configurations. 

One can use the lowest relevant post-Newtonian ex- 
pansion to support (qualitatively) our results. A look at 
the Post- Newtonian equations of motion [35] shows that 
the leading spin-orbit interaction is of 1.5PN order, while 
the leading spin-spin interaction is of 2PN order. Kid- 
der [351 ] gives the radiated energy and momenta of the 
point-particle spin-orbit and spin-spin couplings up to 
2PN order. The relevant equations for our nearly head- 
on collisions are: 

^ 8m 4 2 ft k \ 

E N = — yz— j^y, 15a 
15r 4 

^ = --^{-110^ + 452^}, (15b) 
105r 4 I r J 

E so = 0, (L N = 0), (15c) 

3?7T) 2 ^ 

Ess = ^{Sx ■ S 2 )v 2 r , (15d) 
5r° 

where r is the binary separation, 2m = mi ± v r is 
the relative velocity, N denotes Newtonian terms, PN de- 
notes (non-spin-orbit and non-spin-spin) post-Newtonian 
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FIG. 5: The final remnant horizon mass versus time for 
the HS++ configuration. The top panel shows the horizon 
mass for the three resolutions with grid-spacings h — M/25, 
h = M/30, and h — M/40, and the two Richardson extrapo- 
lations based on leading third-order and leading fourth-order 
errors. The bottom panel shows the differences in the horizon 
mass for the h = M/25 and h = M/30 runs as well as the dif- 
ferences in the mass for the h = M/30 and h — M/40 runs. 
This latter difference is rescaled by 1.57052 to demonstrate 
3.7-order convergence. The accurate extrapolations (as evi- 
dent by the agreement of the two extrapolations) to infinite 
resolution allows for a sufficiently precise determination of the 
radiated energy to model its dependence on the initial spins. 



terms, SO denotes spin-orbit coupling terms, and SS de- 
notes spin-spin coupling terms. Note that the lowest- 
order spin-correction to the radiated energy is quadratic 
in S. Likewise we can obtain the radiated angular mo- 
mentum 

J N = 0, [L N = 0), (16a) 
J PN = 0, (L N = 0), (16b) 
Jss = {a N ~ a S s ~ f), (16c) 

— * 4 777,^" -* — » 

J SO = x {Vr X (Si + S 2 )) 

( 9 2fh 1 4m 2 1 
L r r z J 

which displays the leading linear dependence on the spin, 



FIG. 6: The angular momentum of the final remnant black 
hole for the HS++ configuration measured by the isolated 
horizon method. The top panel shows the horizon spin for 
the three resolutions with grid-spacing h — M /25, h — M/30, 
and h = M/40, as well as the two Richardson extrapolations 
based on leading third-order and leading fourth-order errors. 
The bottom panel shows the differences in the horizon spin 
for the h = M/25 and h = M/30 runs as well as the differ- 
ences in the spin for the h — M/30 and h — M/40 runs. This 
latter difference is rescaled by 1.54 to demonstrate 4th-order 
convergence. The accurate extrapolations (as evident by the 
agreement of the two extrapolations) to infinite resolution al- 
lows for a sufficiently precise determination of the radiated 
angular momentum to model its dependence on the initial 
spins. 



and the radiated linear momentum 

Pn — 0, (Sm = mi — m 2 — 0) (17a) 

P PN = 0, (Sm = mi - m 2 = 0) (17b) 

P ss = (at this order) (17c) 
l, Sm 3 

P so = ~Y^ Vr ^ r X A) = °' 

(A = 2(5i - S 2 ) = 0). (17d) 

Note that equation (|17d|) gives a non-vanishing linear 
momentum radiated for anti-parallel spins. The spin- 
orbit interaction is the only one that reverses signs when 
spins change signs. This explains the different stability 
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FIG. 7: The gravitational energy radiated up to time t as 
measured by an observer at r = 35M for the 1 = 2 modes 
of the HSH — h and MS++ configurations (for our configura- 
tions 7?2 2 = R2-2, and 72 2 = The energy from the 
spurious initial burst leads to the first plateaus centered at 
t = 60M. These initial bursts of energy, while significant, 
are dominated by the total emitted radiation. We multiply 
h±2 from the MS++ configuration (S = 0.1M 2 ) by 4.3 to 
demonstrate that quadratic (plus smaller higher-order terms) 
dependence of 72 ±2 on S. 



and radiation properties of parallel and anti-parallel con- 
figurations. We plan to report the results of this spin up 
- spin down case leading to radiated linear momentum in 
a forthcoming paper. 

These spinning-black-hole-binary scenarios were also 
studied in the complementary regime of the close limit 
approximation by Gleiser and Dommguez [40j ] . They 
examined the behavior of non-boosted, slowly rotating 
holes to second perturbative order. In that approxima- 
tion the second-order terms proportional to the spin and 
the square of the proper distance, I, give rise to radiation 
of the angular momentum, J, proportional to the spin 

J rad w 0.0023 (2S) {l/Mf. (18) 

The analysis also predicts that the radiation of energy 
can be expanded as 

^ « 7.8 x 10- 4 {2S/M 2 ) 4 + 9.8 x 10~ 5 (Z/M) 4 

- 2.8 x 10~ 5 {2S/M 2 ) 2 {l/M) 2 

+ 1.3 x 10~ 3 (2S/M 2 ) 2 (l/M) i , (19) 

The first term, being proportional only to S" 4 , is gen- 
erated by the radiation content of the BY initial data. 
The second term, proportional to I , is the radiation 
for two non-spinning masses in a head-on collision. The 
third term is second-perturbative-order spin-orbit cou- 
pling term, and is the one that leads to the radiation of 
angular momentum. 



V. DISCUSSION OF RESULTS 

We have studied the spin-up of initially non-spinning 
and initially corotating black-hole binaries during the last 
stages of the inspiral orbital motion and found that it is 
much smaller than that needed to achieve tidal locking. 
This means that the time scale of tidal effects to spin up 
the black holes is much longer than the time scale for 
gravitational radiation to drive the merger (and hence 
increase the orbital frequency). This also implies that 
in an astrophysical scenario, black holes with low spin 
merge with essentially the same spin magnitude they had 
when the binary formed (note that Eq. ([7]) implies that 

\S\ = 0). We also note that it is this short timescale for 
gravitational radiation during the last few orbits that is 
responsible for the absence of an ISCO (innermost sta- 
ble circular orbit) when the two holes have comparable 
masses. The concept of the ISCO arises naturally in the 
extreme mass ratio limit, where there is a cutoff between 
stable and unstable orbits. In this limit, the instability 
of orbits inside the ISCO is not driven by gravitational 
radiation, but by non-radiative terms in the particle's 
equation of motion. Thus the particle very slowly inspi- 
rals due to the emission of weak gravitational radiation, 
and then plunges when it reaches the ISCO. On the other 
hand, in the comparable mass regime, there is no such 
cutoff because gravitational radiation dominates the dy- 
namics during the entire late-time inspiral and plunge 
phases. 

In Fig. [5] we show the tracks of the horizon centroids 
in the SO and SC configurations. Although coordinate 
dependent, the tracks help illustrate the lack of an ISCO 
in the comparable mass regime. Rather than showing 
an abrupt change near the last orbit, the tracks show a 
smooth increase in plunge velocity over the entire evolu- 
tion. Note that the SC track is a tighter spiral than the 
SO track, which is consistent with the results found pre- 
viously by the authors [U]. Additional horizon tracks 
for multi-orbit non-spinning binaries can be found in 
Refs. [4, 5], while tracks for highly-spinning binaries can 
be found in Ref. [ll|. 

After the merger, the final remnant black hole formed 
from the SO configuration has a specific spin of a/Mn = 
0.688 ± 0.001 (where Mh is the remnant horizon mass) 
while remnant formed from the SC configuration has a 
specific spin of cl/Mh = 0.717±0.001. Using these values 
in conjunction those found in Ref. [ll[ (see Table IVlj) . we 
find (from a least-squares fit) that the specific spin of the 
final remnant horizon (cl/Mh)\r varies with the initial 
specific spins of the two (equal-mass, equal-spin) black 
holes (a/m)\i according to: 

(a/M H )\ R = 0.6879 + 0.2952 ({a/m)\i) 

- 0.0374 ((a/m)|/) 2 , (20) 

with an error of ±0.001 in the range (—0.757 < 
(a/m)\i < 0.757). From this fit we extrapolate that 
the remnant horizon will have a specific spin no smaller 
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TABLE VI: The measured remnant horizon specific spin 
(a/M_f/)|ij and energy radiated E ra d, as well as the predicted 
remnant horizon specific spin (a/MH)\ P red and energy radi- 
ated E ra d/M\ pre d (based on a least-squares fit) for quasi- 
circular, equal-mass, equal-spin binaries with initial specific 
spins (a/m)\i. 



(a/m)|j (a/A/ H ) 



(a/Mn )\ P red E rad /M E rad /M\ pr 



-0.757 0.443 ±0.001 0.4430 

0.000 0.688 ±0.001 0.6878 

0.1001 0.717 ±0.001 0.7169 

0.757 0.890 ±0.001 0.8900 



(2.2 ±. 01)% 2.20% 
(3.5 ±0.1)% 3.48% 
(3.8 ±0.1)% 3.79% 
(6.7 ±0.2)% 6.70% 



-1.0 
±1.0 



0.355 
0.946 



2.2% 
8.1% 



-4 1 1 1 1 1 1 1 1 1 

-4 -2 2 4 

x/M 

FIG. 8: The horizon tracks for the SO and SC configurations 
with tick marks every 10M of evolution (we display only one 
of the SC tracks for clarity). Note that the SO and SC tracks 
agree well for the first 100M, and that after 100M the SC 
tracks show a tighter spiral and orbit longer before merging. 
The SO configuration completes slightly less than 1.75 orbits, 
while the SC configuration completes slightly more. 



than 0.355 (for the case of initially-maximal, anti-aligned 
spins) and no larger than cl/Mh — 0.946 (for the case of 
initially-maximal, aligned spins), which supports our as- 
sertion in (llj | that the cosmic censorship conjecture is 
not violated by the merger of two black holes. A least- 
squares fit of the energy radiated versus the initial indi- 
vidual spins yields 

= 0.0348 ± 0.0297 ((a/m)|/) ± 

0.0170 ((a/m)|/) 2 , (21) 

with an error of ±0.001 in the range (—0.757 < 
(a/m)\i < 0.757). Based on this fit we extrapolate that 
a quasi-circular binary will radiated no more than 8.1% 
of its mass during the final few orbits and merger. 

We have also found that the tidal coupling in black- 
hole binaries responsible for the transfers of angular mo- 
mentum from the individual spins to orbital is relatively 
small. Nevertheless, we have been able to compute this 
transfer with sufficient accuracy to model its dependence 
on the initial spin. We find that the radiated angular mo- 
mentum, which is a measure of the angular momentum 
transferred from spin to orbit, has a linear-plus-cubic de- 
pendence on the initial spins; in qualitative agreement 
both with the post-Newtonian (see Eq. (|16d[) ) and close- 
limit (see Eq. HI])) predictions. The radiated energy has 
a quadratic-plus-quartic dependence on the initial spin; 
again, in qualitative agreement with the post-Newtonian 
(see Eqs. (|15p ) and close- limit (see Eq. (fH?)) ) predictions. 

Unfortunately the speculations of Price and Whe- 



lan [4l[ suggesting this spin-up effect might be strong, 
and even dominant, are not confirmed by our computa- 
tions. We see neither the possibility of a binary stalling 
by this effect nor the dominance of the I = 3 mode sug- 
gested in Ref. (III ]. 



APPENDIX A: ISOLATED HORIZON SPIN AND 
HORIZON CIRCUMFERENCES 

The most common technique for evaluating the spin of 
a numerically evolved black hole is to use the ratio of the 
polar to equatorial circumferences (42J. For a Kerr hole 
this ratio C r — C v /C eq is given by 

c 1 + y^Rgg gf (°/m) 2 ) 

{ (l + v /l-(a/m) 2 )V' 
(Al) 

where a/m (m being the horizon mass) is the specific 
spin and E(x) is the complete elliptic integral of the 
second kind. After a common horizon forms, C r dis- 
plays quasinormal ringing with a constant offset. To 
obtain the specific spin in finite time one uses a non- 
linear least-squares fit of C T versus time to C r = C'o ± 
C 1 exp[~C 2 t/M}sm[C 3 t/M + C 4 ] to obtain C , • • • , C 4 
(Co being the asymptotic value of C r ), and then a/m is 
obtained by inverting Eq. (|Aip . The authors have used 
this technique to obtain highl y-a ccurate evaluations of 
the remnant spin in Refs. [^Tdlllll]. However, in the pre- 
merger case this technique is not accurate. There are 
two reasons for this. First, one cannot separate the large 
oscillations in C r from the secular growth because the 
behavior of the oscillations is not know a-priori. Thus, 
there is a large error in the calculated spin. Second, the 
distortions on the horizon that lead to C r ^ 1 are due 
to both the spin as well as tidal distortion of the hori- 
zon, and when the black holes are close, the tidal distor- 
tion dominate. Thus even if one can obtain an accurate 
measurement of the secular part of C r , one would then 
need to disentangle the spin component of C r from the 
tidal distortion component. Finally, the horizon circum- 
ference calculations are not coordinate invariant and lead 
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FIG. 9: The specific spin of the individual horizons for the SC 
orbiting binary configuration as calculated using the horizon 
circumferences and the isolated horizon technique. Here a xz 



is obtained from C xz /C eq via Eq. (|A1|I . 



from C yz /C e , 



Note the large oscillations and secular growth in a xz and a yz 
compared to the much more accurate isolated horizon spin. 



to errors when the coordinate are distorted (i.e. when the 
coordinate circles used in the calculations do not corre- 
spond to the correct coordinate circles in Boyer-Lindquist 
coordinates). 

In Fig. [5] we show a/m for the SC orbiting-binary con- 
figuration calculated using the isolated horizon method 
and the circumference method. In the plot we show a/m 



calculated using the ratios of the xz and yz polar cir- 
cumferences to the equatorial (i.e. xy) circumference, as 
well a/m derived from the isolated horizon technique. 
Note that the initial values for the spin 0.03 and 0.07 for 
the two circumference-based calculations are significantly 
smaller than the expected 0.1, while the isolated horizon 
technique gives the correct value. The difference between 
the spin calculated from the two circumferences is a mea- 
sure of the error. In this case that error is larger than 
the effect we want to measure. Note also that the spin 
derived from the circumferences shows a much stronger 
increase with time. This strong increase is driven by the 
increasing tidal distortion of the horizon (as can be in- 
ferred from the increase in amplitude of the oscillations). 
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